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We imaged reflectance and variable fluorescence in 25 cyanobacterial mats from 
four distant sites around the globe to assess, at different scales of resolution, 
spatial variabilities in the physiological parameters characterizing their photosynthetic 
capacity, including the absorptivity by chlorophyll a (-Achi), maximum quantum yield of 
photosynthesis (Vmax), and light acclimation irradiance (4). Generally, these parameters 
significantly varied within individual mats on a sub-millimeter scale, with about 2-fold 
higher variability in the vertical than in the horizontal direction. The average vertical profiles 
of Vmax and decreased with depth in the mat, while /\chi exhibited a sub-surface 
maximum. The within-mat variability was comparable to, but often larger than, the 
between-sites variability, whereas the within-site variabilities (i.e., between samples from 
the same site) were generally lowest. When compared based on averaged values of 
their photosynthetic parameters, mats clustered according to their site of origin. Similar 
clustering was found when the community composition of the mats' cyanobacterial layers 
were compared by automated ribosomal intergenic spacer analysis (ARISA), indicating 
a significant link between the microbial community composition and function. Although 
this link is likely the result of community adaptation to the prevailing site-specific 
environmental conditions, our present data is insufficient to identify the main factors 
determining these patterns. Nevertheless, this study demonstrates that the spatial 
variability in the photosynthetic capacity and light acclimation of benthic phototrophic 
microbial communities is at least as large on a sub-millimeter scale as it is on a global scale, 
and suggests that this pattern of variability scaling is similar for the microbial community 
composition. 
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INTRODUCTION 

A major challenge in microbial ecology is to understand how 
the structure, composition, and function of microbial commu- 
nities are linked, how microbial communities are influenced by 
environmental conditions, and how they contribute to the local 
and global cycling of elements. Examples of microbial com- 
munities that have been extensively studied from all of these 
perspectives are cyanobacterial mats, which are highly com- 
pacted microbial ecosystems consisting of diverse phototrophic 
and heterotrophic populations (van Gemerden, 1993; Stal, 2000). 
The interest in their study stems from the generally accepted 
assumption that they represent a modern analog of the earli- 
est complex ecosystems on Earth (Seckbach and Oren, 2010). 
Additionally, their vertically stratified and compact structure 
offers good possibilities for studying microbial and biogeochemi- 
cal interactions in well-controlled laboratory mesocosms. 



Environments that harbor cyanobacterial mats vary greatly 
with respect to parameters such as nutrient concentrations, 
temperature, salinity, or input of mineral particles. This vari- 
ability leads to differences in the structure of the mats with 
respect to the density and distribution of microbial cells, con- 
tent of inorganic particles, and characteristics of exopolymers 
that bind the mat matrix together. In spite of these differ- 
ences, one feature that these environments have in common is 
their extremity with respect to at least one of the environmen- 
tal parameters, typically salinity, temperature, or pH (Caumette 
et al., 1994; Garcia-Pichel et al, 1999; Des Marais, 2003; Abed 
et al., 2006; Seckbach and Oren, 2010). This environmental con- 
straint is essential for excluding, or at least minimizing, the 
influence of grazing and bioturbation, which would otherwise 
disturb the perennial growth and laminated structure of the 
mats. 
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A typical feature of cyanobacterial mats are steep vertical gra- 
dients of physical and chemical parameters such as light, O2, pH, 
and H2S (van Gemerden, 1993; Stal, 2000). While the gradients 
in light intensity and spectral composition are the consequence 
of strong and wavelength-dependent absorption and scattering 
by photosynthetically active (photopigments) and inactive (min- 
eral particles, organic detritus) components in the mat (Kiihl and 
Jorgensen, 1992, 1994; Kiihl et al., 1994), the steep chemical gra- 
dients form due to mass transfer limitation in a volume that 
is densely packed with active microbial cells (Kiihl et al., 1996; 
Wieland and Kiihl, 2000a; Jonkers et al., 2003; Garcia de Lomas 
etal, 2005). 

A critical process in cyanobacterial mats is the photosynthetic 
activity of the cyanobacterial population, which occurs in the 
uppermost layer of the mat (so-called euphotic zone) and sup- 
ports diverse heterotrophic populations in the mat ecosystem 
through the production of organic substrates and O2 (Ntibel 
et al., 1999; Roeselers et al, 2007). Cyanobacterial photosynthe- 
sis depends on a number of environmental parameters, including 
temperature, salinity, concentration of nutrients, intensity and 
spectral quality of light, and exposure to H2S (Kiihl, 1993, 2005; 
Kiihl and Fenchel, 2000; Wieland et al, 2003; Pinckney et al, 
2011). Since these parameters vary between different sites that 
harbor cyanobacterial mats as well as within the mats them- 
selves (see above), it is expected that the rate and efficiency of 
photosynthesis in the mats will exhibit strong geographical as 
well as micrometer-scale variability. Although previous measure- 
ments demonstrated that photosynthetic rates and efficiency are 
strongly variable within mats (Kiihl et al., 1994; Al-Najjar et al., 
2010), presently it is not known how this within-mat variability 
compares to the variability between mats from different envi- 
ronments. Additionally, it is not known whether there is a link 
between the photosynthetic capacity in the euphotic zone of the 
mats and the composition of the corresponding microbial com- 
munity, and to which extent are these properties of the mats 
determined by the parameters characterizing their environment. 

To address these issues, we compared 25 samples of cyanobac- 
terial mats collected from four distant geographical locations 
(United Arab Emirates, Australia, Brazil, and Spain) with respect 
to their photosynthetic capacity and microbial community com- 
position. Our focus was on the cyanobacterial layer at the top 
of the mats, which was assumed to be a good proxy for the 
mats' photosynthetically active zone. We used absorptivity by 
chlorophyll a (Achi), maximum quantum yield of photosystem 
11 (iniax), and light acclimation irradiance (Ik) as parameters 
characterizing the photosynthetic capacity and adaptation of the 
cyanobacterial populations. These parameters were measured 
with a sub-millimeter spatial resolution across vertical sections 
of the cyanobacterial layers using hyper-spectral and variable 
chlorophyll fluorescence imaging. Differences in the microbial 
community composition of the cyanobacterial layers were quanti- 
fied by automated ribosomal intergenic spacer analysis (ARISA). 
Possible links between the parameters characterizing the pho- 
tosynthetic capacity, microbial community composition, and 
environmental parameters were identified using multivariate sta- 
tistical methods. We hypothesized that the microbial community 
composition and photosynthetic capacity in the cyanobacterial 



layers are linked, and that the photosynthetic potential and light 
acclimation of the cyanobacterial populations are more strongly 
influenced by the steep vertical gradients within the mat ecosys- 
tem than by the environmental parameters characterizing their 
habitat. 

MATERIALS AND METHODS 
SAMPLES 

The studied cyanobacterial mats originated from four sites: an 
intertidal flat near Abu-Dhabi, UAE (AD mats), an intertidal flat 
in the Exmouth Gulf in Australia (AU mats; Lovelock et al., 2010; 
Adame et al., 2012), the hypersaline lake Lagoa Vermelha in Brazil 
(BR mats ; Vasconcelos et al., 2006), and the hypersaline lake La 
Salada de Chiprana in Spain (SP mats; Jonkers et al., 2003). The 
mat samples were collected between 2003 and 2008 in at least two 
replicates from each site and incubated under artificial illumina- 
tion (lOh light/ 14 h dark cycles, wavelength range 400-700 nm) 
and at approximately constant temperature until the measure- 
ments, which were conducted in 2009. More details about the 
collection sites and incubation conditions are given in Table 1. 
During the incubation, the appearance of the AD, AU, and BR 
mats did not change, whereas the SP mats gradually changed from 
thinly laminated structures (see Jonkers et al., 2003) to thicker 
structures featuring a loosen upper layer composed of a mixture 
of exopolymers and suspended particles. 

MEASUREMENT PROTOCOL 

First, the mat samples were pre-incubated for about 12 h at room 
temperature and incident irradiance of 100|xmol photons roT^ 
s^^. Subsequently, they were vertically sectioned and immedi- 
ately afterwards variable chlorophyll fluorescence and spectral 
reflectance were measured to characterize, respectively, the photo- 
synthetic potential and pigments in the mats. This was done with 
a high spatial resolution (~20 [im) across the vertical sections of 
the mats using imaging cameras (see below). Immediately after 
imaging, cyanobacteria-dominated layers close to the mats sur- 
face were cut-off with a sterile scalpel and prepared for the ARISA 
(see below). Identification of these layers was based on their char- 
acteristic dark-green appearance, microscopic observations and 
the results of the imaging analyses, which were rapidly obtained 
by the image processing routines developed during this study (see 
below). 

PAM IMAGING OF THE VARIABLE CHLOROPHYLL FLUORESCENCE 

Pulse amplitude modulated (PAM) imaging of the variable 
chlorophyll fluorescence was done with the Imaging-PAM system 
(Walz GmbH, Germany), using red light- emitting diodes for the 
excitation of the chlorophyll a fluorescence from cyanobacteria. 
For each mat, a vertical section of the mat sample was placed 
on its side in a Petri dish and covered with a few millimeters 
of seawater (salinity of 32, temperature 15°C). After 15min of 
dark adaptation, images of the minimum (_Fo) and maximum 
(Frn) fluorescence yields in the dark-adapted state were recorded. 
Subsequently, rapid light curves (RLC) (Schreiber et al., 1996, 
1997) were measured by increasing the actinic irradiance from 
0 to 1700 |xmol photons m^^ s^' in time intervals of 3 min for 
each irradiance and acquiring the fluorescence yield images under 
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Table 1 | Characteristics of the sampling sites, collection details, and incubation conditions for the studied microbial mats. 



Site abbreviation 



AD 



AU 



SP 



BR 



ENVIRONMENTAL PARAMETERS" 



Type of site 
Salinity range (%o) 
Temperature (°C) 
Inundation period 
Water depth when 
inundated (m) 
Type of substrate 
N (mg/l) 
P (mg/l) 

SAMPLE COLLECTION 

Year 

Geographic location 
Replicates 

Within-site replicates 
separation (m) 

INCUBATION CONDITIONS 

Salinity range (%o) 
TrC): 

Downwelling irradiance'' 



Intertidal flat 
35-200 
25-55 
Few hours everyday 
1 

Very fine grained sediment 
2 
0.03 

2007 

24°31'20"N, 
54°26'50"E 
13 
-50 



For 



Intertidal flat 

35-50 

24-28 
~6 h per month 
2 



Very fine grained sediment 
840 
96 

2008 
22M3'7"S, 



Hypersaline lake 
78-90 
15-28 
Continuous 
1.5 

Very fine grained sediment 
1.2 
0.04 



Lagoon (mixed with seawater) 
40-150 
25-29 
Continuous 
1-2 

Organic material 
NA" 
0.23 



114°34'6"E 



2003 

41°14'30"N, 
0°10'50"W 
2 



35-200 
26-28 

380= 



35-50 
26-28 
120-180' 



80-90 

20 
240S 



{[imo\ photons m 



1): 



2008 
22°53'15.2"S, 
42°6'38.14"W 
2 

-10 



35-200 
26-28 
380== 



" Information extracted from literature: Ai-Najjar et al. (2010) for AD mats; Lovelock et al. (2010) for AU mats; Jonkers et al. (2003) for SP mats: Vasconcelos et al. 
(2006) for BH mats. 
^ Value not measured. 

Taken from an aquarium where the mats were grown for~8 years. 
'^Quantified by a calibrated PAR quantum irradiance sensor (LI-190 Quantum) connected to a light meter (LI-250, both from LI-COR Biosciences). 
" Light source: AQUALINE 10000, MH 400W, Germany 
^ Light source: Envirolite, UK. 

^ Light source: cool white fluorescent tubes T8 (32W), Philips, Germany. 



actinic illumination {F') and during the saturating pulse (F^') at 
the end of each interval. The intensity and duration of the satu- 
rating pulse was ~2400 (imol photons m^^ s^^ and 0.8 s, respec- 
tively, and the irradiance levels of the actinic light were calibrated 
using a PAR (400-700 nm) quantum irradiance sensor (LI-190 
Quantum) connected to a light meter (LI-250, both from LI-COR 
Biosciences) positioned in the same place as the mat sample. As 
the sample was relatively small and it was lying on its side, irradi- 
ance was evenly distributed across the vertical section of the mat. 

IMAGING OF CHLOROPHYLL A ABSORPTIVITY 

In addition to variable fluorescence imaging, the Imaging-PAM 
system was used to image reflectance of the mats in the red (J^r) 
and near-infrared (Rnii) region. These images were used to cal- 
culate the chlorophyll a absorptivity as A^hi = (^nir — Ri)IRnk^ 
which was taken as a proxy for chlorophyll a concentration in 
the mats. Because the spectral resolution of these measurements 
was insufficient, chlorophyll a absorptivity was additionally quan- 
tified by hyper-spectral imaging (Kiihl and Polerecky, 2008; 
Polerecky et al, 2009). Specifically, each mat sample in the Petri 
dish was placed on a motorized stage, illuminated with a halogen 
bulb (Philips, type 6423) emitting in the visible to near-infrared 



range (400-900 nm), and scanned with a hyper-spectral imaging 
system (VNIR-100, Themis, Themis Vision, USA). Spectral nor- 
malization was achieved by scanning a gray reference standard 
with 40% reflectance (SRS-40-020, Labsphere Inc., USA). After 
verification that the reflectance spectrum had a characteristic chl 
a edge in the wavelength range of 700-720 nm and a flat plateau 
above 720 nm (see, e.g., Polerecky et al, 2009), the image of chl a 
absorptivity was calculated as A^hi = (Riso — ^675)/J^750) where 
i?750 and i?675 are the reflectance images measured at 750 nm 
(no absorption by chl a) and 675 nm (maximal chl a absorp- 
tion), respectively. No difference was found between the chl a 
absorptivity values obtained by hyper-spectral imaging and by the 
Imaging-PAM system (data not shown). Therefore, the latter were 
used in the subsequent analysis, as they allowed perfect alignment 
with the variable fluorescence images. 

DNA EXTRACTION AND ARISA 

DNA from the cyanobacterial layers in the studied mats was 
extracted and purified using the UltraClean soil DNA isolation 
kit (MO BIO Laboratories, Inc., Carlsbad, CA, USA) according 
to the manufacturer's instructions. For each mat sample, PCR 
(50 (il) were conducted in triplicates and contained 1 x PCR 
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buffer (Promega, Madison, WI, USA), 2.5 mM MgCl2 (Promega), 
0.25 mM of 40 mM dNTP mix (Promega), bovine serum albumin 
(3 (ig/|xl, final concentration), 25 ng extracted DNA, 400 nM each 
of universal primer ITSF (5'-GTCGTAACAAGGTAGCCGTA- 
3') and eubacterial ITSReub (5'-GCCAAGGCATCCACC-3'; 
Cardinale et al., 2004) labeled with the phosphoramidite dye 
HEX, and 0.05 units GoTaq polymerase (Promega). All sub- 
sequent steps, including the PGR protocol, purity of the PGR 
products, labeling of the products, discrimination of the PCR- 
amplified fragments via capillary electrophoresis (ABI PRISM 
3130xZ Genetic Analyzer, Applied Biosystems) and the subsequent 
statistical analysis of ARISA profiles (quality control, binning, 
merging) were done as previously described (Boer et al., 2009; 
Ramette, 2009). 

PROCESSING AND ANALYSIS OF THE FLUORESCENCE YIELD IMAGES 

Using the fluorescence yield images, the quantum yield of PSII in 
the dark-adapted state was calculated as Yq = (^^m — Po)/Fm and 
the effective quantum yield of PSII at a given actinic irradiance, 
7 > 0, as 7 = (fm' — P')IFui'- Subsequently, the values of Y were 
plotted as a function of I to determine the maximum effective 
quantum yield of PSII (denoted as Ymax) and the actinic irradi- 
ance at which Y reached Ymax (denoted as /max)- Based on the 
theoretical background of the saturation pulse method (Baker, 
2008), Ymax represents the maximal photosynthetic potential of 
the cyanobacterial population in the mat. In plants, Y decreases 
monotonously with I and this maximal potential is reached in 
the dark-adapted state, i.e., Ymax = io and Jmax = 0 (White and 
Critchley, 1999). However, in our measurements with cyanobac- 
terial mats the Y-I relationship was not monotonous and these 
parameters were typically related as Ymax > and /max > 0 
(Figure 1), which is why we additionally determined also /max- 

In photosynthesis research, light acclimation is determined 
from the relationship between the rate of photosynthesis, P, and 
irradiance, 1. This so-called P-I curve is close to linear at low irra- 
diance levels and approaches saturation at high irradiance levels. 
The actual value of the light acclimation irradiance (denoted here 
as 7k) depends on the mathematical model that describes it, and 



is typically obtained by fitting the P-I data (see, e.g., Piatt and 
Jassby, 1976). To obtain Ik from our data, we assumed that the 
parameter Y measured by the saturation pulse method is given by 
Y = a P/I, where a is a proportionality constant whose value is 
not important in this study (but see Campbell et al., 1998), and 
considered the following model to describe the P-I relationship: 

PU) = Pmax[l - exp(- I/Ik)] - 8Y/a7exp(- I/I^). (1) 

The first term in this model is equivalent to that proposed 
by Webb et al. (1974), whereas the second term was necessary 
to describe the non-monotonous relationship between Y and / 
observed in this study (Figure 1). Using Equation (1) and assum- 
ing that 7m ^ ^k) which was the case in our measurements, Y can 
be written approximately as 

Y ^ Ym(/k//)[l - exp(- //7k)] - 8 Y exp(- ///max), (2) 

where Ym = Ymax/(1 - W/4)> 8Y = Ym - Yq and /max ~ 

if <SC (see Figure IB). Thus, by determining the values of 
Yq, Ymax, and /max and fitting the rest of the measured Y-I val- 
ues with the model in Equation (2), it was possible to determine 
the remaining fitting parameter /k. Data processing required for 
this analysis, including quantification of the images of Ymax, 
■fmax, and /k, their variabilities across the vertical and horizon- 
tal directions as well as average vertical and horizontal profiles, 
was done in Matlab (The Math Works Inc., Natick, MA) using 
the program Look@PAM developed during this study. This pro- 
gram is available on the internet (http://www.microsen-wiki.net/ 
pamimaging:lookatpam) . 

STATISTICAL ANALYSES 

The significance of univariate response data as a function of 
categorical factors was tested using one-way analysis of vari- 
ance (ANOVA), after verifying the normality (Shapiro-Wilk nor- 
mality test) of the response variable at p = 0.05. Community 
differences were visualized by non-metric multidimensional scal- 
ing (MDS) ordination based on the Bray-Curtis dissimilarity 




FIGURE 1 I (A) Examples of quantum yields of photosystem II, X in 
the studied cyanobacterial mats, as measured by pulse amplitude 
modulated (PAM) imaging at different irradiances, /. Symbols and error 
bars represent, respectively, the mean and SD calculated from 5x5 



pixels in the image. (B) Mathematical model describing the V vs. / 
relationship. The main features of the relationship are annotated. The 
corresponding relationship between photosynthesis, P, and / is also 
shown. 
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matrix between samples, and significance of community differ- 
ences between groups of samples was determined by Analysis 
of Similarity (ANOSIM) tests. Prior to analyzing ARISA profiles 
conjointly with functional or environmental variables, a consen- 
sus community profile was obtained for each sample by merging 
the triplicate ARISA PGR and by considering an OTU present if it 
appeared at least twice among the triplicates (Ramette, 2009). The 
merged table was Hellinger-transformed to minimize the effects 
of the strongly right-skewed distribution curve (Ramette, 2007). 
To assess the link between the microbial community composition 
and function, the Procrustes superimposition approach was used 
to estimate the concordance of scores originating from two inde- 
pendent ordinations after rotating, translating, and dilating one 
of them, while keeping the other ordination coordinates constant 
(Gower, 1975). Significance of the rotation statistic was assessed 
by Monte-Garlo permutations (Peres-Neto and Jackson, 2001). 
All statistical tests were carried out with the statistical platform _R 
(http://cran.r-project.org/) and multivariate community analyses 
with the vegan package (Oksanen et al., 2012). 

RESULTS 

AH studied mats had a clear laminated structure, with a charac- 
teristic dark-green layer at or close to the mat surface in each of 
them (see examples of true color images in Figure 2). The AD 
and AU mats had a brown layer and a red layer underneath the 
dark-green layer, while the deepest layer was black. Additionally, 
the AD mats were covered by a thin orange gelatinous layer. The 
BR mats showed a different layering pattern with an upper dark- 
green layer followed by a thicker dark-pink gelatinous layer. The 



SP mats had a similar structure as the AD mats, except the distinct 
layers were thicker and the surface gelatinous layer had whitish to 
light-brown appearance. 

Hyperspectral imaging revealed that the dark-green layer in 
the mats had a pronounced absorption at wavelengths corre- 
sponding to the maximal absorption by chlorophyll a (675 nm) 
and phycocyanin (625 nm; Figure SI). Since these pigments are 
characteristic for cyanobacteria, we could be confident that the 
dark-green layer contained abundant cyanobacterial populations 
and could therefore be referred to as the cyanobacterial layer. This 
conclusion was supported by microscopic observations (data not 
shown). 

VARIABILITY OF PHYSIOLOGICAL PARAMETERS ON DIFFERENT 
SPATIAL SCALES 

Imaging of A^hi, Vmax! Vax) and 7k in mats collected from differ- 
ent sites made it possible to investigate how the variability in these 
physiological parameters changes depending on the scale at which 
it is determined, including the micrometer-scale (within-mat 
variability), meter-scale (within-site variability), and global-scale 
(between-site variability). The within-mat variability was quanti- 
fied in three ways: as standard deviation (SD) of the values from 
the entire image of the cyanobacterial layer and as SD of the aver- 
age vertical and horizontal profiles. The within-site variability was 
calculated as SD of the average values over the cyanobacterial layer 
in each mat from the respective site, while the global variability as 
SD of the average values for each site. 

As demonstrated by the images, all studied physiological 
parameters exhibited pronounced micrometer-scale variability 
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within the cyanobacterial layer of the studied mats (Figure 2). 
Two-factorial analysis of variance performed on individual 
images of A^hi, ymaxi and Ik. using vertical and horizontal posi- 
tion in the mat as factors, revealed that the percentage of the total 
within-mat variance explained by the vertical position (25-30% 
for Achi, 45-50% for Y^^^, and 30-35% for 4) was about 1.5 to 
2-fold larger than the percentage of variance explained by the hor- 
izontal position (15-20% for A^hi, 15-30% for y^ax, 15-25% for 
/k). This was consistent with the comparison of the SD for the 
average vertical and horizontal profiles, which showed that the 
former were about 2-fold larger than the latter (Figures 3A-C). 
Thus, on average, the vertical variability in Achi, Ymax> and 7k was 
about twice as high as the horizontal one. A significant portion of 
the within-mat variance was explained by the interaction between 
the vertical and horizontal position (55-60% for A^u, 30-35% 
for Ymax! and 45-50% for 7k), consistent with the clearly visi- 
ble variation of the vertical profiles of these parameters along the 
horizontal direction for each individual mat sample (Figure 2). 
With respect to Imax^ the vertical and horizontal variabilities 
were similar (Figure 3D), each explaining about 10% of the total 
within-mat variability, while the remaining 80% was explained 
by their interaction. The within-mat variability of 7mi,x in the AD 
mats was significantly larger than in the other mats (Figure 3D), 
which was mainly because the I^m values were generally larger in 
the AD mats. 

Average vertical profiles of A^hi for AD, BR, and SP mats were 
characterized by a subsurface maximum, whereas a monotonous 
decrease with depth was observed for AU mats (Figure S2A). 
Average Ymax and 7k generally decreased with depth (Figures 
S2B,C) and were significantly correlated (p < 0.001) for each mat 



sample. Spatial trends in 7inax were not so clear, having patchy dis- 
tributions for some mats while on average increasing with depth 
for others (Figure 2 and Figure S2D). Depending on the mats, 
vertical variation across the cyanobacterial layer, expressed as a 
coefficient of variation (SD/mean) of the average vertical profile, 
reached 10-30% for Achi, 5-45% for y^ax. 10-60% for 7k, and 
5-100% for 7inax (Figure 3). 

The within-site variability of A^hi calculated for the mats from 
the AD and AU sites was lower than the within-mat variabilities 
calculated for the individual mats from these sites (Figure 3A). 
This was often the case also for Ymax. ^k) and 7max; but some- 
times the within-mat variability was lower than the within-site 
one (Figures 3B-D). Such comparison could not be done reli- 
ably for the BR and SP mats because of the limited number of 
replicates {N = 2) for these two sites. 

The global-scale variability was larger than the within-site 
variability for all measured physiological parameters (Figure 3). 
In contrast, in many cases it was lower than the within- 
mat variability. This relationship was most pronounced for 
Achi and 7k, where it was observed in 21 (for A^hi) and 
12 (for 7k) out of 24 mat samples (Figures 3A,C). Clearly, 
this was primarily due to the pronounced variability in the 
vertical direction, whereas the horizontal within-mat variabil- 
ity was almost always lower than the global-scale variabil- 
ity. On the other hand, the global-scale variability in Ymax 
and 7max was mostly larger than the within-mat variability, 
except for a few AD mats (Figures 3B,D). The global-scale 
variability, expressed as a coefficient of variation, was about 
15% for Achi, 40% for Ymax, 30% for 7k, and 80% for 7max 
(Figure 3). 
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FIGURE 3 I Variabilities of the physiological parameters Ach\ (A), Vn 
(B), /k (C), and /^ax (D) calculated at different spatial scales, 
including the micrometer (within-mat), meter (within-site), and 



global (between sites) scales. See text for details. The variabilities are 
expressed as standard deviations, SD (left axis), and as a coefficient of 
variation (right axis). 
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CLUSTERING OF MATS BASED ON AVERAGE PHYSIOLOGICAL 
PARAMETERS 

Physiological parameters AcU, Ymax) Vax; and Ik. when averaged 
over the cyanobacterial layers, varied significantly between the 
sampled sites (Figure 4, Table 2). For example, AU and BR mats 
had, on average, about 2-fold larger Ymax than AD and SP mats, 
SP mats had the lowest 7k> whereas AD mats had the largest values 
of /max- When a distance matrix was calculated from the aver- 
age values of Achi, Ymaxj and /max using Euclidean metric, its 
visualization in an MDS plot revealed clear clustering of the mats 
according to the site of their origin (Figure 4C). 

CLUSTERING OF MATS BASED ON THE MICROBIAL COMMUNITY 
COMPOSITION 

ARISA fingerprinting revealed that the microbial communities in 
the mats from a given site were more similar to each other than 
to those from other sites (Figure 5A). This marked endemism 
was further supported by significant ANOSIM test (P = 0.0001). 
Overall, the sampled bacterial communities shared between 18 
and 34% OTUs, with 84 out of 398 OTUs (21.1%) found every- 
where (i.e., at least in one mat from a given site) and 314 OTUs 
being mat-specific. 

LINK BETWEEN COMMUNITY COMPOSITION AND FUNCTION 

To assess the degree of concordance between community com- 
position and their potential function (photosynthetic capacity), 
we compared the bacterial community composition as depicted 
by the NMDS plot in Figure 5A with that of the configuration 
of Achi, ymax> ^k. and /max (Figure 4C) by Procrustes analyses. 
A significant concordance between ordinations was found (r = 
0.605, P = 0.002, based on 1000 permutations; Figure 5B), sug- 
gesting that distinct communities were associated with distinct 
(potential) functions. 

ENVIRONMENTAL PARAMETERS 

Most parameters characterizing the habitat of the studied mats 
were obtained from literature and are summarized in Table 1 . The 



sites strikingly differed with respect to nutrient concentrations in 
the overlying water, which were extremely high for the AU mats. 
This was very likely because the sediments at this site contained 
large biomass of insect larvae (alive and dead), which could be 
a source of nutrients (Behie et al., 2012). Additionally, the site 
was inundated in brief and infrequent intervals (during a spring 
tide). Another notable difference between the sites was due to the 
salinity fluctuations, which were very large at the AD and BR sites. 

DISCUSSION 

SCALING OF VARIABILITY IN PHYSIOLOGICAL PARAMETERS 

The imaging approach used in this study enabled us to assess how 
the spatial variability in parameters characterizing the photosyn- 
thetic potential in cyanobacterial mats scales depending on the 
sampling resolution (from sub-miUimeter to thousands of km). 
Our main result is that this scaling pattern differed depending on 
the studied parameter, with chlorophyll a absorptivity (Achi) and 
acclimation intensity (Ik) having on average the largest variability 
on the micrometer-scale (mostly in the vertical direction) while 
the maximum quantum yield of PSII (Ymax) and the irradiance 
at which this maximum yield is reached (/max) were most vari- 
able on the global scale (between sites). This contrast suggests that 
the different physiological parameters are controlled by different 
environmental factors, as discussed below. 

The pronounced micrometer-scale vertical variability in Achi 
was most likely due to the combined effects of light and nutrients. 
The growth of photosynthetic microbial mats is typically limited 
by nutrients, with most nutrients assimilated in the cyanobacte- 
rial layer originating from organic matter remineralization under- 
neath the layer (lonkers et al, 2003). On the other hand, light 
quantity strongly attenuates with depth due to intense absorp- 
tion by photopigments and abiotic components of the mat matrix 
(Kuhl and Jorgensen, 1992, 1994; Kuhl et al, 1994). Thus, because 
of the opposing gradients in light and nutrients, cyanobacterial 
growth at the top and bottom of the cyanobacterial layer is likely 
limited by nutrients and light, respectively, whereas a location 
with an optimal supply of light and nutrients is somewhere in the 
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FIGURE 4 I (A,B) Relationships between physiological parameters A^^], /max, 
/k, and /max in the studied cyanobacterial mats. Symbols represent averages over 
the cyanobacterial layers, thin and thick error-bars depict the total and vertical 



variability (expressed as standard deviation, SD) within the layers, respectively. 
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Table 2 | Average values and standard deviations of the physiological parameters characterizing the mats from the different sites. 





AD 


AU 


BR 


SP 






0.66 ±0.06 [2]*= 


0.80±0.08[1] 


0.72 ±0.05 [1-2] 


0.57 ±0.10 [2] 


0.0004 


Knax 


0.13±0.05 [2] 


0.25±0.03[1] 


0.280±0.002[1] 


0.13±0.01 [2] 


8 X 10-6 


/k(|xmol ) 


85±23[1] 


88±16[1] 


66 ±11 [1-2] 


32 ±5 [2] 


0.01 


/max(|J.mol m-2 s-^ 


24±12[1] 


5 ±3 [2] 


2.4 ±0.2 [2] 


6±4[2] 


1 X 10-5 



^Probability that the means between the different sites are equal, as determined byANOVA. The /max values were log-transformed before ANOVA to ensure variance 
homogeneity. 

^ Means ranked as 111 are significantly larger than means ranked as [2], means ranked as 11-2] are not significantly different from those ranked as [1] and 12]. 



Procrustes errors 




FIGURE 5 I (A) Variation in tine microbial community composition within 
the cyanobacterial layers of the studied mats. A Bray-Curtis dissimilarity 
matrix was calculated based on ARISA community profiles and is 
displayed in a 2D ordination space (associated stress value of 0.182). 
Grouping lines were added for each site a posteriori to highlight the site 
specificity in the community patterns. (B) Procrustes analysis of the link 



Dimension 1 



between the composition and function of the microbial communities in 
the studied mats. Open symbols represent the ARISA-based NMDS data, 
while the closed ones represent the functional data (Figure 4). Note the 
small rotation angle between the two sets of ordination axes. AD mats 
(m2-ml5); BR mats (m17-m18); AU mats (m19-m30); SP mats (m35 
and m38). 



middle. This suggests that the cyanobacterial biomass should have 
a maximum somewhere around this optimal location. Assuming 
that the measured chlorophyll a absorptivity, Achi, is a proxy for 
cyanobacterial biomass, our results are consistent with this inter- 
pretation: while Achi had a clear subsurface maximum in mats 
from sites where the N and P concentrations in the overlying 
water were low (AD, SP, BR), it decreased sharply from the surface 
in the mats from the AU site, where the nutrient concentrations in 
the overlying water were very high (Table 1, Figure S2). This latter 
characteristic of the AU site was most likely also the main factor 
responsible for on average the highest values of Achi (Figure 4A, 
Table 2) and its vertical variability (Figure 3A) in the AU mats 
as compared to the mats from the other sites. Additional explana- 
tion for the observed within-mat variability in Achi is the tendency 
of phototrophic cells to have a larger cellular pigment content 
when grown at lower irradiances (Falkowski, 1980; Kirk, 2011). 
Thus, at least in the mats from AD, SP, and BR, the observed 
subsurface maximum in Achi could additionally be due to this 
adaptation of cyanobacterial cells to light limitation that progres- 
sively increases with depth in the mat. Together, this suggests that 
the average cyanobacterial biomass in the cyanobacterial layer of 



the studied mats was mostly determined by the overlying water 
nutrient content, whereas the micrometer-scale distribution of 
the biomass, and possibly also of the average pigment content in 
the cyanobacterial cells, was additionally shaped by light. 

With respect to the light acclimation intensity, light avail- 
ability was likely the most important factor that determined both 
its micrometer- scale and global-scale variation. Assuming that 
cyanobacteria adapt to local light conditions, the sharp decrease 
in light intensity with depth in mats should result in cyanobac- 
terial populations in deeper parts of the euphotic zone being 
acclimated to lower light intensities. This is what we generally 
observed for mats from all sites (Figure S2). Although we do not 
have reliable data on site-specific downwelling irradiances, the 
observation that the global-scale variability of Ik was comparable 
to, and in many occasions lower than, the micrometer-scale vari- 
ability suggested that light variation within the photosynthetically 
active cyanobacterial layer was similar or larger than the variation 
in the average light dose received by the mats from the different 
sites. 

Although the micrometer-scale variability of the maximum 
effective quantum yield, Yniaxi was strongly positively correlated 
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with the Hght acdimation intensity for each studied mat (see e.g., 
Figure S2), light availability is not a likely factor that controls its 
spatial distribution. This is primarily because Ymax is, by defini- 
tion, a proxy for the maximum photosynthetic potential, which 
is reached at light intensities that are considerably lower than 
the acclimation intensity (in the dark for plants, at /max <K for 
the cyanobacterial populations in this study; Figure 1). Since the 
quantum yield of PSII determined by the PAM measurement of 
variable chlorophyll fluorescence relates to the redox state of the 
plastoquinone (PQ) pool in the photosynthetic electron transport 
chain (Allen, 2003), parameters that affect this redox state are 
likely affecting also Ymax- In the context of cyanobacterial mats, 
H2S and O2 are likely candidates. 

The response of cyanobacterial photosynthesis to H2S is well- 
documented (Cohen et al, 1986; Jorgensen et al., 1986; Miller 
and Bebout, 2004), with several species being able to perform 
an oxygenic photosynthesis using H2S as the electron donor. In 
such cyanobacteria, H2S oxidation is facilitated by the activ- 
ity of sulfide-quinone-reductase (SQR) (Bronstein et al., 2000; 
Griesbeck et al., 2000), which reduces the PQ molecule using 
the electrons transferred from H2S and can therefore decrease 
the apparent quantum yield of PSII. Our preliminary experi- 
ments with an axenic cyanobacterial culture embedded in agarose 
showed that exposure to H2S in the |xM range lead to a rapid 
decrease in the quantum yield Y, both in the dark and light, and 
the yield recovery in the light occurred only after H2S decreased 
below a threshold in the [xM range (Figure S3). A similar negative 
effect of H2S on Y was observed in other cyanobacterial systems 
(e.g., living stromatolites; Kromkamp et al, 2007). Therefore, the 
pronounced decrease in Ymax with depth in the cyanobacterial 
layer may be due to the steep increase in H2S in this layer at low- 
light conditions, which is a well-documented phenomenon in 
mats (e.g., Wieland and Kiihl, 2000b; Jonkers et al., 2003; Garcia 
de Lomas et al, 2005; Al-Thani et al, 2014). 

The effect of O2 on Y is likely because photosynthetic and 
respiratory electron transport chains in a cyanobacterial cell 
share three electron carriers in the thylakoid membrane (PQ, 
cytochrome bef and the terminal oxidases; Vermaas, 2001), which 
makes O2 an important electron acceptor for electrons from the 
reduced PQ (chlororespiration) and also for those coming from 
PSI via ferredoxin (Mehler reaction) (Schreiber et al., 2002). 
Presence of O2 can therefore contribute to a partial oxidation of 
the PQ pool, which could manifest itself as an increase in the 
apparent Y. The steep decrease in O2 with depth, which is typi- 
cal for cyanobacterial mats under low light conditions (Wieland 
and Kuhl, 2000a, 2006; Al-Najjar et al, 2012), could therefore 
additionally be responsible for the observed decline of Ymax with 
depth. 

In addition to environmental parameters mentioned above, 
part of the observed spatial variability of Ymax> 4 and A^hi 
could be due to artifacts linked to autofluorescence (for Ymax 
and 7k) and absorption (for A^hi) by photosynthetically inactive 
components in the cyanobacterial layer, such as dead/inactive 
cyanobacterial cells, pigment degradation products or mineral 
particles. Such components are abundant in the cyanobacte- 
rial layer of microbial mats (e.g., Kiihl et al., 1994; Al-Najjar 
et al., 2012), and because their fluorescence is not variable but 



possibly comparable to that from the PSII of photosyntheti- 
cally active cells, the detected quantum yield of PSII can appear 
lower. Similarly, their absorption in the same wavelength range 
as chlorophyll a could increase the apparent A^hi. These effects 
may have contributed to the differences in the "baseline" of these 
parameters and thus also their apparent global variability. The 
importance of this contribution could, however, not be estimated 
based on our present data. 

In contrast to the monotonous decrease in the effective quan- 
tum yield, Y, with the actinic irradiance, I, which is charac- 
teristic for plants (White and Critchley, 1999) and eukaryotic 
algae (Flameling and Kromkamp, 1998), our measurements in 
cyanobacterial mats showed that the relationship between Y and I 
was not monotonous (Figure 1). Instead, Y reached a maximum 
at actinic intensities /max that varied in the range from one to 
several tens of |xmol photons m^^ s^^ (Figure 4B). The spatial 
patterns for /max within a given mat were not very clear (see, e.g.. 
Figure 2), primarily because of the noise in the fluorescence yield 
images. Although we were able to find a suitable mathematical 
expression for this relationship (Equation 2), further research is 
required to understand why Y reached a maximum at intensities 
/max > 0 and what controls this behavior. 

IMPLICATIONS FOR MICROBIAL ECOLOGY 

Our imaging data demonstrate that the variability of a spe- 
cific function (here photosynthesis) in a microbial community 
(here the cyanobacterial layer) can be at least as large on the 
micrometer- scale within the community as it is on the global 
scale between communities from different locations. As discussed 
above, this is primarily due to the strong within-mat variabil- 
ity of the physico-chemical parameters that control the func- 
tion. Furthermore, our comparison of the cyanobacterial layers 
based on their average functional (photosynthetic) potential and 
microbial community composition revealed that the mats clus- 
tered according to their site of origin (Figures 4C, 5A) and that 
the functional and compositional data were significantly linked 
(Figure 5B). Based on these results we expect that if physico- 
chemical parameters exhibit pronounced variability within a 
microbial community, which is common in transport-limited 
systems such microbial mats and sediments, the community com- 
position will also vary at least as much on the micrometer-scale 
(within community) as it does on the global scale (between com- 
munities). This implies that analyses of the microbial community 
composition that ignore this micrometer-scale variability (e.g., by 
sampling bulk volumes of sediments rather than layer-by-layer) 
may have a limited value in identifying correlations between the 
microbial composition, function and the corresponding envi- 
ronmental settings. In other words, when sampling microbial 
communities in order to identify these correlations, the prefix 
"micro" should refer not only to the size of the inhabitants but 
also to the required spatial resolution of sampling. 

Interestingly, despite the large micrometer-scale variability of 
the functional potential within each cyanobacterial layer com- 
munity, distinct communities from a given site were on aver- 
age characterized by a distinct, albeit highly variable, functional 
potential (Figure 4). This suggests that out of the environmental 
parameters discussed in the previous section, there are possibly 
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one or more factors that have the most dominant influence on 
the composition and function of the studied cyanobacterial mats. 
Unfortunately, these dominant factors could not be identified 
from the limited dataset presented in this study. 
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